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For interdependent networks with identity dependency map, percolation is exactly the same with that on a 
single network and follows a second-order phase transition, while for random dependency, percolation follows a 
first-order phase transition. In real networks, the dependency relations between networks are neither identical 
nor completely random. Thus in this paper, we study the influence of randomness for dependency maps on 
the robustness of interdependent lattice networks. We introduce approximate entropy(Api?n) as the measure 
of randomness of the dependency maps. We find that there is critical ApEric below which the percolation is 
continuous, but for larger ApEn, it is a first-order transition. With the increment of ApEn, the Pc increases 
until ApEn reaching ApEn[, and then remains almost constant. The time scale of the system shows rich 
properties as ApEn increases. Our results uncover that randomness is one of the important factors that lead 
to cascading failures of spatially interdependent networks. 
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The interdependent networks which fully con¬ 
sider the interactions between networks have 
been used to model real complex systems better. 
Robustness is one of most important properties 
for interdependent networks especially spatially 
interdependent networks, since most of the in¬ 
frastructure networks are spatially networks. In 
real interdependent networks, the dependency re¬ 
lationship are not usually random. Thus we an¬ 
alyze how the randomness of dependency map 
affects the robustness of interdependent lattices 
which are used to model the spatially interdepen¬ 
dent networks. We find that the randomness of 
dependency map between networks is quite criti¬ 
cal for the robustness of interdependent lattices. 


I. INTRODUCTION 

Robustness is one of the most important properties of 
complex networks and has been widely explored on sin¬ 
gle networks in the last decadei^— . However, complex 
systems are rarely isolated. The more casual situation is 
that networks usually interact with other networks such 
as transportation networks and hnancial systems^rJ^. hr 
the case of interdependent networks, conclusions are of¬ 
ten far different from single networks. In particular, a 
removal of a very small fraction of nodes can lead to 
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catastrophic failures on the whole networbi^. A theo¬ 
retical framework based on percolation theory has been 
established to analyze the resilience of interdependent 
systems2,ii^, and much details have been exploredi^^— . 
The fraction of interdependent nodes is one important 
factor that will influence the phase transition of the 
networks2ii^. Also, the overlap of links can significantly 
change the properties of the percolation, and there is a 
critical point above which the emergence of the mutu¬ 
ally connected component is continuous^. The presence 
of degree correlations in multiplex networks can modify 
drastically the percolation threshol d^^'^^ . 

Most previous models have focused on interdependent 
random and scale-free networks in which space restric¬ 
tions are not considered. However, many real-world sys¬ 
tems such as power grid networks and computer networks 
are embedded in two-dimensional spac o^^'^^ . In interde¬ 
pendent random and scale-free networks, the overlap of 
links and degree correlations will change the properties 
of phase transition. Nevertheless for spatially embedded 
interdependent networks which are modeled as square 
lattices, the overlap of links or the degree correlations of 
nodes lose their significance, since their network topolo¬ 
gies are identical. The spatially interdependent networks 
are extremely vulnerable. Any fraction of interdepen¬ 
dent nodes will lead to first order transitions^. From an 
identical dependency map to totally random dependency 
map, the randomness of the dependency map may be one 
of the most important factors leading to the emergence 
of discontinuous percolation. In most real interdepen¬ 
dent systems, dependencies are neither totally random 
nor identical. Research on the resilience of intermedi¬ 
ate systems that he somewhere in between these two ex- 
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tremes is of high practical significance and needs further 
exploration. 

From this perspective, we study the relationship be¬ 
tween the dependency’s randomness and stability of the 
system of two interdependent spatially embedded net¬ 
works. We use approximate entropy(>lpi?n) as the mea¬ 
sure of randomness. One of the big challenges here is how 
to introduce controlled degree of randomness into the sys¬ 
tem. Therefore, we propose an intermediate model which 
describes the system with dependency map between iden¬ 
tical map and totally random map. Inspired by the 
constructing procedure of the Watt-Strogatz small-world 
model^, starting from an identical dependency map, we 
rewire each dependency link at random with probability 
q. By increasing q from 0 to 1, the ApEn increases mono- 
tonically. Therefore, the traverse of randomness can be 
generally represented by q. We reveal that there is a 
critical value qc, for which the percolation transition be¬ 
comes continuous, whereas for any q > qc the collapse 
is discontinuous. Changing the topologies on a single 
layer, we discover that qc is different for interdependent 
scale-free networks, Watts-Strogatz networks, and Erdos- 
Renyi networks. The critical threshold increases with q 
when q < q'c and remains approximately constant when 
q > q'c- Additionally, we present an analytical method 
for time scale of cascade failures based on critical p and 
find that the four topologies display rich transient prop¬ 
erties when q changes from 0 to 1. Finally, we analyze 
the influence of limited dependency length on spatial net¬ 
works. We show that with the same dependency length, 
a linearly dependent system is always continuous, but not 
continuous for some locally randomly dependent system. 
Our results show that the randomness of dependency may 
be one of important factors for extreme vulnerability of 
spatially interdependent systems. 


II. MODEL DESCRIPTION 

Our model of interdependent networks is realized via 
two networks(A^ = 10®) A and B under full dependency. 
Here one network is the copy of the opposite network 
and their average degree < k >= 4(the same as a square 
lattice). The degree distribution of the scale-free network 
is < fc where A = 3.0. In each network, each node 
has two types of links: connectivity link and dependency 
link. Also every node in network A is connected with one 
and only one node in network B. For a square lattice, each 
node is connected to its four nearest neighbors within the 
same lattice via connectivity links. All dependencies in 
our model are mutual and bidirectional. Dependency is 
taken to mean that if a node in network A is removed 
from the system and a node in B that depends on it 
will be removed from B as well. Thus failures of nodes 
iterate until mutually connected giant component of both 
networks emerges. This process is called cascade failures 
and see Methods for details of cascade process of the 
system. 






FIG. 1. The interdependent square lattice with the rewiring 
probability of dependency links q = 0 , 0.25, 0.50, 1.00 re¬ 
spectively. When q = 0, the dependency map is an identical 
mapping i.e. node i in network A is dependent on node j 
in network B where i = j(EigT.a) . When q = 1.00, the 
situation is the same as totally random mapping(Fig.l.b). 
When q = 0.25 or 0.5 , the situation is between both ex- 
tremes(Fig.l.c,d). 


There are two extreme situations, i) node f in A de¬ 
pends on node j in the B such that j = i. We call it 
identity dependency map(FiglT]a). ii) The random de¬ 
pendency map as most papers considered(Fig.l.b). Like 
the constructing procedure of the Watt-Strogatz small- 
world model, starting from the identity dependency map, 
we rewire each dependency link at random with proba¬ 
bility q, while guaranteeing that each node in A depends 
on one and only one node in B(0 < g < 1). We sample 
q = 0, 0.25, 0.50,1 and plot them in Fig.l. 

Note: We must hgure out that our model is differ¬ 
ent from partially interdependent lattices proposed by 
Amir Bashan et al2^. In partially interdependent lat¬ 
tices, there are interdependent lattices with a fraction q 
of interdependent nodes and the remaining 1 — g of nodes 
autonomous. In our model, however, the remaining 1 — g 
nodes are connected with the identical nodes in opposite 
network. The subtle difference between these two models 
will lead to huge difference in cascade failure process. We 
will illustrate it in Figl^ In FigI2 we can see that the 
cascade failures process differs much between these two 
models: with the same g = 5/9 and p = 4/9, the size of 
the giant component in our model is 0/9, while the size of 
giant component in partially interdependent networks is 
4/9. This apparently show that the our model is different 
partially interdependent lattices. 


III. RESULTS 

Entropy can be used to measure the randomness 
effectively^^. In fact, approximate entTopy(ApEn) is 
adopted in this paper for computation convenience. 
When g = 0, ApEn is nearly 0, and when g = 1, it reaches 
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FIG. 2. Difference of cascade failures between partially randomly interdependent lattices and partially interdependent lattices. 
The blue points stand for the survived nodes, while the red points stand for the failure nodes. Figure sequence a-f stands for 
the cascade failures process in partially randomly interdependent lattices with q — 5/9(fraction of nodes that are randomly 
dependent and the remaining 1 —q of nodes are dependent with the identical nodes in the opposite network) and p = 4/9(fraction 
of nodes initially removed). Figure sequence A-C stands for the cascade failures process in partially interdependent networks 
with q = 5/9(fraction of nodes that are dependent and the remaining 1 — q are autonomous)and p = 4/9. It can be obviously 
seen with the same q, the size of giant component in partially randomly interdependent lattices is 0/9, while the size of giant 
component in partially interdependent lattices is 4/9. The cascade failures process differs for these two models. 


its maximum. The ApEn is a continuous and monoton- 
ically increasing function of (7(Fig.3). However the ran¬ 
domness is not fully represented by rewiring dependency 
links, since the locally randomly interdependent lattice^l 
in which q = 1 is not totally randomly interdependent 
but with length constraints. Then considering a more 
casual situation, the permutation of 1 ^ fV cannot be 
exhausted by rewiring the dependency links of identical 
mapping at probability q. But as the approximate en¬ 
tropy changes continuously with q, we can change q to 
generally represent all approximate entropies. 

Through simulation, we find that there is a critical 
qc ~ 0.13 for an interdependent lattice network be¬ 
low which the percolation is second-order but hrst-order 
above. In Fig.4, we can see that for q = 0.1, the de¬ 
crease of giant component occurs in multiple steps, char¬ 
acteristic of a second-order transition. For q = 0.2 and 
q = 1.0, the giant component may completely collapse 
by removal of a small fraction of nodes, characteristic 
of a first-order transition(Fig.4). pc increases with q. 
When q is relatively small, pc increases approximately 
linearly. And when q > q'^, pc remains almost con¬ 
stant, i.e. as more random the dependency is, as more 
fragile is the system. Analogously, for interdependent 
scale-free networks, Watts-Strogatz networks and Erdds- 
Renyi networks, there is also a critical qc- We find that 
qf < qf^ < q^^ < qf^. Additionally pc of lattice net¬ 
works is generally greater than that of other networks. 
Smaller fraction 1—pc will lead to cascade failures. This 
means that an interdependent scale-free network is most 
robust under random attacks, while the square lattice is 
most vulnerable. A random network is more stable than 



FIG. 3. The value of ApEn under different q. When q 
= 0, ApEn is nearly 0,and when q = 1, ApEn reaches its 
maximum. The ApEn is a continuous function of q and it 
change monotonously as the increment of q. 


a small-world one. 

The time scale of cascade failures, i.e. the time that the 
interdependent networks need to collapse to the station¬ 
ary state is one evidently important merit for system’s 
resilience. The number of iterations( NOI) increases and 
reaches its peak at the critical qc and goes quickly down 
to a small value with p (Fig.5). So the NOI at pc is 
an effective measure for time scale of the system. The 
NOI at Pc (be.Pc) is a function of g. p® increases quickly 
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FIG. 4. Relations of the size of Poo at steady state after ran¬ 
dom failure of a fraction 1 — p of the nodes on interdependent 
square lattices. The green circles, red squares and blue tri¬ 
angles for q — 0.1,g = 0.2 and q = 1.0 respectively. The 
numerical results are obtained by averaging 100 realizations 
of networks consisting oi N = 10®. For q = 0.1, the phase 
transition of the system is continuous since the giant com¬ 
ponent occurs at multiple steps. For q — 0.2 and q — 1.0, 
however, the transition is first-order as the giant component 
collapses even by removing a very small fraction of nodes. 



FIG. 5. The function of number of iterations(NOI) vs p in 
interdependent lattice network(g = 1). The numerical results 
are obtained by averaging 100 realizations of networks con¬ 
sisting of = 10“* nodes. The vertical red line is the critical 
line. On the left side of it, the system collapses down(blue 
circle), while a giant component remains functional on the 
right side(green circle). There is a sharp divergence of the 
NOI when p approaches Pc- 


with q when q < qc and declines very gently above qc 
for interdependent lattice networks. For a scale-free net¬ 
work, it increases until qc and then starts to decline. For 
random networks and small-world networks, it increases 
monotonously with q but the variation tendency becomes 
gentle above qc(Fig.6). All four interdependent systems 
have variation tendency’s changes around qc- pi of inter¬ 
dependent lattice networks is greater than those of scale- 
free, small-world and random networks when q < q^}. pi 
of interdependent square lattice is smaller than those of 
all other three network types when q > q^^. The NOI re¬ 
flects the time scale of system collapse. Our results show 
that the transient characteristics of the four systems go 
through rich changes with the variation of q. 

Finally we check locally interdependent networks in 
which the distance between two interdependent nodes is 
limited(d < r i.e. \xi — X 2 \ < r and j/i — ?/ 2 | < r in 
reference^) . For simplicity, we consider one more special 
condition. Here we split the whole network into r=i=r small 
blocks and dependency links are randomized within each 
block. We find that there is critical distance Tc ~ 25 un¬ 
der which the percolation is continuous but discontinuous 
above rc(Fig.7). The Tc is greater than in^i because 
the randomness(approximate entropy) here is lower than 
that in^Iwith the same distance. The corresponding ap¬ 
proximate entropy ApEn ~ 0.923 of is approximately 
equal to ApEuc- Compared with locally random depen¬ 
dency, the linear dependency map is more robust. For the 
linear dependency map, the percolation is always contin- 
uous(Fig.7). Although the dependency distance changes 


strongly, the approximate entropies of those dependency 
maps are almost equal to 0. So their percolation prop¬ 
erties are nearly the same as that of a single lattice. It 
is thus clear that the randomness may be a more impor¬ 
tant factor leading to cascade failures than dependency 
distance. 

Furthermore, it is possible that the randomness of de¬ 
pendency is related to other metrics of interdependent 
networks such as dimension. The dimension of net¬ 
works is a function of the distribution of link lengths^^. 
For spatially embedded networks, the dimension is one 
of the most fundamental quantities to characterize its 
structure and very likely will influence its percolation 
property^. However, how interdependency relationships 
between networks change the dimension of the system 
has not been figured out. In reference^!, the authors dis¬ 
cover that the dependency length play an critical role 
in the percolation transition. However, we found that 
under linear dependency map, the change of dependency 
length influence the percolation property little. From the 
discrepancy of those two situations, it can be inferred 
that local property of dependency relationship make a 
big difference. And the local property of dependency will 
directly influence the local topological inter-similarity be¬ 
tween networks. Randomness happen to reflects the lo¬ 
cal property of dependency (we can see this from the 
computation steps of approximate entropy in Methods). 
In spatially interdependent embedded networks, the lo¬ 
cal characteristics of dependency can be more intuitively 
characterized as the relative length of dependency links. 
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FIG. 6. The change of NOI at pc with q. The brown, cyan, red 
and blue circles stand for NOI at pc of interdependent square 
lattice, scale-free network, Watt-Strogatz network and Erdos- 
Renyi network respectively. For interdependent random net¬ 
works and small-world networks, NOI increases monotonously 
with q. However for interdependent lattice network and scale- 
free network, there is one critical q^ ■ When q < q^ , the 
NOI increases approximatively linearly, and when q > qc , 
the NOI starts to decreases. 


Under linear dependency map, the relative length of de¬ 
pendency links and the approximate entropy of depen¬ 
dency map are nearly 0. No matter how large r is, they 
changes little and are still nearly 0, so the percolation 
changes little. On the other hand, the smaller the rela¬ 
tive length is , the less dimension is changed from that of 
single network. There should exit some relations between 
dimension and the randomness of the dependency map. 


IV. DISCUSSION 

In many real interdependent systems such as coupled 
power grid and communication network, the dependency 
relationship is neither completely regular nor completely 
random, but lies somewhere between these two extremes. 
The transition from regular to random dependency is one 
of the keys to extreme vulnerability of spatially interde¬ 
pendent systems. From the proposed intermediate cas¬ 
cade failure model (from regular to random dependency), 
we find there is a transform from continuous percolation 
to discontinuous percolation with the randomness vari¬ 
ation of the dependency map between two interdepen¬ 
dent networks. We emphasize the generic character of 
our model because the dependency map could influence 
not only the resilience but also synchronization, disease 
spreading and other dynamic process in interdependent 
networks. With suitable modification, our model could 
be applied to understand the dynamical process in most 
real interdependent systems since the dependency maps 



FIG. 7. The fraction of nodes in the giant component as a 
function of q. For locally random interdependent network, 
when r = 8, the system represents the characteristic of a 
second-order transition. For r = 25 and r = 50, the gi¬ 
ant component may completely collapse by removal of even 
a single additional node, which represents the characteris¬ 
tic of a first-order transition. However for linearly interde¬ 
pendent networks, the transitions are second-order one when 
r = 8, 25, 50 and even r = L. 

between networks are more various and complicated in¬ 
stead of totally random dependency or regular depen¬ 
dency. 

The time scale of cascade failures is evidently impor¬ 
tant for system’s resilience, but it has received little at¬ 
tentions in the analysis of resilience so far. In different 
dynamic processes, the characteristic time scales of sys¬ 
tems vary greatly. For instance biological systems, social 
systems and financial market dynamics have time scale 
much longer than that of cascade failures of power grid. 
Our analytic method based on critical p is simple and 
effective for characterizing the time scale of different sys¬ 
tems. Generally, the system which has a shorter time 
scale demands higher requirements for our responding 
speed to catastrophe and brings much bigger challenges 
for us to take mitigation actions than those with longer 
time scale. Therefore, our method may provide a clue for 
research on revealing the transient mechanism and miti¬ 
gation of cascade failures in interdependent networks. 


V. METHODS 
A. Approximate entropy 

The randomness of the dependency maps the interde¬ 
pendent square lattice is measured by approximate en¬ 
tropy. Entropy can effectively reflects the randomness 
of a sequence. However, for computation convenience, 
we choose the approximate entropy as the measure of 
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the randomness of the system. The approximate en¬ 
tropy is denoted by ApEn and is computed by following 
steps^Sii^^: 

A. Given a series 

X{i) = [u{i), u{i -I- 1), • • • , u{i -l-m — l)],i = l^Af — 
TO -I- 1 

B. Count the distance between the vector X{i) and 
other X{j) for each i: 

d[X{i), X(j)] = max \u{i + k) — u{j + k)\ 

k—O^m — 1 

C. Given an threshold, count the the ratio between 
the number such that d[X{i),X{j)] < r for each i and 
the number of the vectors ie.N — m+ 1 (devoted by ("(r). 
as: 

_ {t/te number of(d\X{i),X(j)])<r} 

'-'r V ) N-m+l 

Generally, C^{r) reflects m-dimensional pattern 

N — 771+1 

D+’-(r) = E InCrir) 

E. ApEn{m,r) = $"*(r) — <i)™+i(r) 

Parameter selection: here we choose to = 2 , and r = 
0.2 * (standard deviation of u). 


B. Percolation transition 


The percolation transition is studied by randomly re¬ 
moving a fraction 1 — p of nodes and the links attached to 
them from both networks simultaneously. Then, on each 
network, clusters which are detached from the largest 
connected component are removed. After that, the nodes 
in each network which lost their supporting nodes in 
the opposite network are removed. This in turn causes 
more clusters to break off from the giant component and 
this process is continued until no more clusters break 
off. First, we analyze the situation with totally ran¬ 
dom dependency maps. After the initial attack, only 
Pi = Poo{p) remains functional. Each node in A that is 
removed causes the removal of its interdependent node in 
B. Then only Poo(pi) nodes in B remain alive. This pro¬ 
duces further damage in A and causes cascading failures. 
The cascade failures can be represented by the recursive 
equation for the survived fraction 


Po =P, 

p 

Pi = ^Poo(po) = Poo{p), 
VO 

Pi = -^Poo(Pi-l)- 

-M-1 


( 1 ) 


In the limit i —>■ oo, Eq.(l) yields the equation for the 
mutual giant component at the steady state. 


X = ^/pPoo(x). 


( 2 ) 


Equation (2) can be solved graphically by the intersec¬ 
tion between the curve y = pPoo(x) and the straight 
line y = X. Next, we consider the mutual percolation 
for more casual situations where the dependency is not 


totally random. For every dependency link, there is a 
probability q to rewire it at random. This is equivalent 
to the situation with a fraction 1 — q nodes mapping to 
itself and the remaining q nodes have a random depen¬ 
dency map. The case oiq = 1 corresponds to the scenario 
of a random dependency map and g = 0 is identical to 
the conventional percolation on a single lattice. For the 
initial attack which destroys a fraction 1 — p of nodes, 
[to = (1 — p)A^J nodes are removed. We compute the 
probability Psame that one node in A depends on the 
same node in B. For n nodes in totally random depen¬ 
dency networks, the number of nodes E(n) in the same 
location of both networks is^: 

71 

^ m* CffD(n — to) 

E(n) = ^-, (3) 

n\ 

71 

D(n) = *nl/k\. (4) 

fc =2 

When n is very large, the computation of D(n) is very 
inconvenient. For computation simplication, when n >2, 
we have: 


D{n) ^ [n\/e + 0.5\. (5) 

where e is the Euler’s number and [xj is the integer part 
of X. Then, 


E{n) = 


E (;+^L(n-TO)!/e + 0.5)J 


nl 


771=1 


nl 


( 6 ) 


= E 


1 


1 


—'^e(TO—1)! 2(n — to)!(to — 1)! ’ 

1 




771 = 0 


2{n — m + 1)!to! 2 


So for each node, the probability that it is in the same 
location of A and B is: 

Psame = (I - q) * P + P * . (7) 

q* N 

When n ^ oo , Psame —^ (1 ~ 9 ) * P- The initial attack 
causes some number of nodes to be disconnected from 
the giant component in both networks A and B. Further¬ 
more, because of the dependency relationship, the nodes 
disconnected from A will lead to further damages. Poo 
increases with Psame- The greater Psame is, the more 
nodes disconnected from the giant component of A over¬ 
lap the nodes in B. So further damage decreases and cas¬ 
cade failures are weakened (or prevented) from the be¬ 
ginning . For g = 0, the cascade failures are prevented 
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from the start and the percolation is continuons. When 
q — I, the totally random dependency map will lead to 
a first-order transition. As q increases, the percolation 
transition changes from a continuous transition to a dis¬ 
continuous one. There must be one critical qc beyond 
which the percolation transition become discontinuous. 
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